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Consider a critical nearest neighbor branching random walk on 
the d-dimensional integer lattice initiated by a single particle at the 
origin. Let G„ be the event that the branching random walk sur- 
vives to generation n. We obtain limit theorems conditional on the 
event Gn for a variety of occupation statistics: (1) Let V„ be the max- 
imal number of particles at a single site at time n. If the offspring 
distribution has finite ath moment for some integer a > 2, then in 
dimensions 3 and higher, Vn = Op{n^^"); and if the offspring dis- 
tribution has an exponentially decaying tail, then Vn = Op (log n) in 
dimensions 3 and higher, and Vn = Op ( (log n)^) in dimension 2. Fur- 
thermore, if the offspring distribution is non-degenerate then P{Vn > 
iSlogn|G„) — >■ 1 for some 5 > 0. (2) Let Mn{j) be the number of 
multiplicity-^ sites in the nth generation, that is, sites occupied by 
exactly j particles. In dimensions 3 and higher, the random variables 
Mn{j)/n converge jointly to multiples of an exponential random vari- 
able. (3) In dimension 2, the number of particles at a "typical" site 
(that is, at the location of a randomly chosen particle of the nth 
generation) is of order Op (log n) , and the number of occupied sites 
is Op (n/ log n). We also show that in dimension 2 there is particle 
clustering around a typical site. 



1. Introduction. A nearest neighbor branching random walk is a discrete- 
time particle system on the integer lattice Z'^ that evolves according to the 
following rule: At each time n = 0, 1, 2, ... , every particle generates a ran- 
dom number of offspring, with offspring distribution Q = {Qi}i>o', each of 
these then moves to a site randomly chosen from among the 2d + 1 sites at 
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distance < 1 from the location of the parent. We shah consider only the 
case where the branching random walk is critical, that is, where the mean 
number of offspring per particle is 1, and we shall assume throughout that 
the offspring distribution has finite, positive variance cr^. 

By a well-known theorem of Kolmogorov (see Athreya and Ney (1972), 
ch.l) if the branching process is initiated by a single particle, and if G„ is 
the event that the process survives to generation n, then 

(1.1) 7r„ := P(G„) ~ 

Therefore, if the branching random walk is started with n particles at time 0, 
then the number of initial particles whose families survive to time n follows, 
approximately for large n, a Poisson distribution with mean 2/(T^, and the 
number of particles Zn alive at time n is of order Op{n). In fact, in this case, 
under suitable hypotheses on the initial distribution of particles, the measure- 
valued process associated with the branching random walk converges, after 
rescaling, to the super-Brownian motion Xt with variance parameter cr^ (see 
e.g., Etheridge (2000)). In dimensions 2 and higher, the random measure Xf 
is, for each t > 0, almost surely singular with respect to the Lebesgue measure 
on M*^; and when d > 3, the measure Xf spreads its mass over the support in 
a fairly uniform manner (Perkins (1988)), and in fact can be recovered from 
its support (Perkins (1989)). It is natural to conjecture that this uniformity 
also holds, in a suitable sense, for critical branching random walk, and that 
the maximal number of particles at a single site at time n does not grow 
rapidly in n. Our main results show that this is indeed the case. For ease of 
exposition, we will state our results as conditional limit theorems given the 
event G„ of survival to generation n. Corresponding unconditional results 
for branching random walks started by n particles could easily be deduced. 

We shall assume throughout the paper, unless otherwise specified, that 
the branching random walk is initiated by a single particle located at the 

^Allowing particles to remain at the same locations as their parents with positive 
probability eliminates some annoying periodicity problems that would require tedious, 
but routine, arguments to circumvent. Our main results could be proved under much less 
restrictive hypotheses on the jump distribution. 
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origin at time 0. Set 



(1.2) 





= set of particles in generation n; 




= \Zn\ = number of particles in generation n; 


Unix) 


= number of particles at site x in generation n; 




= number of occupied sites in generation n; 


Mn{j) 


= number of multiplicity-j sites in generation n; and 


Vn 


= max Un{x). 



(A multiplicity-j site is a site with exactly j particles.) 

Definition 1. Let X„ be a sequence of random variables, f{n) a se- 
quence of positive real numbers, and Hn a sequence of events of positive prob- 
ability. Say that X„ = Op{f{n)) given Hn if the conditional distributions 
of Xn/f{n) given Hn are tight. Similarly, say that Xn = op{f{n)) given Hn 
if the conditional distributions of X„//(n) given Hn converge weakly to the 
point mass at 0. 

Theorem 2. Assume that the offspring distribution Q has finite ath 
moment for some integer a >2, and that d > 3. Then conditional on Gn, 

(1.3) Vn = Op{n'/''). 

In particular, if Q has finite moments of all orders, then Vn = Op{n^) for all 
e > 0. 

Theorem 3. Assume that the offspring distribution Q has an exponen- 
tially decaying tail, that is, there exists 5 > such that exp(5/) < oo. 
Then conditional on Gn, 

(1.4) K = Op(logn), ifd>3; 

(1.5) K = Op((logn)2), tfd = 2. 

In fact (see Corollary 16 below) for sufficiently large C > the condi- 
tional probabilities P{Vn > Clogn|G„) in dimensions d > 3 and P{Vn > 
C(log n)^|G„) in dimension d = 2 decay polynomially in n. For one-dimensional 
branching random walk, it is known that Vn is of order ^/n (Theorem 7.10 
in Revesz (1994)); stronger results are proved in Lalley (2009). 

Theorem 4. Assume that d >2. Then there exists 5 > 0, depending on 
the offspring distribution Q, such that 

(1.6) lim P{Vn>6logn\Gn) = 1. 

imsart-aop ver. 2009/08/13 file: BRWpublish. tex date: April 8, 2010 



4 



STEVEN P. LALLEY AND XINGHUA ZHENG 



Theorem 3 and Theorem 4 imply that, in dimensions 3 and higher, if 
the offspring distribution has an exponentially decaying tail then Vn is of 
order logn on the event G„ of survival to generation n. In particular, the 
(conditional) distributions of V^/logn are tight, and any weak limit has 
support contained in [(5i,(52] for some 81,82 > (cf. Corollary 16). This 
partly settles an open question (Question 2, p. 79) raised in Revesz (1996). 

Theorem 5. Assume that d > 3. Then conditional on the event Gn, the 
joint distribution of the occupation statistics Mn{j)/n converges as n ^ 00. 
In particular, for certain constants kj such that Yl^iJ ' '^j ~ ^> 

where Y is exponentially distributed with mean 2/a'^. 

This extends the classical theorem of Yaglom, according to which the 
conditional distribution of Zn/n, given that the branching process survives 
to generation n, converges to the exponential law with mean 2/a'^. See 
Athreya and Ney (1972), ch. 1 for a discussion of Yaglom's theorem and 
related results; and Geiger (2000) for an interesting probabilistic proof. 

Theorem 5 implies that in dimensions 3 and higher, most occupied sites 
are occupied by only 0(1) particles. Ultimately, this is a consequence of the 
transience of random walk in dimensions d > 3. Since random walk in di- 
mension d = 2 is recurrent, different behavior should be expected for the 
occupation statistics of branching random walk. In the following theorem 
and throughout this article, we shall use the term typical particle to mean a 
particle chosen randomly from the nth generation Zn. of the branching pro- 
cess (with the choice made independently of the evolution of the branching 
random walk up to time n, according to the uniform distribution on Zn). By 
a typical site we mean the location of a typical particle. 

Theorem 6. In dimension d = 2, the number Tn of particles at a typical 
site at time n is, conditional on the event Gn, of order Op(logn). Moreover, 
for some sufficiently small e > there exists 5 > such that 

(1.8) liminf P(r„ > elognl G„) > 5. 

n— >oo 

We conjecture that the conditional distributions of T„/ log Ti given Gfi con- 
verge in distribution as n — )■ 00. Fleischman Fleischman (1978) has used the 
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method of moments to establish a related result for the number of particles 
at a fixed site at distance 0(1) from the origin. Unfortunately, calculation of 
higher moments for the number of particles at a typical site appears to be 
considerably more difficult, and so the method of moments does not seem to 
be a feasible approach to the conjecture. 

By Yaglom's theorem, conditional on the event of survival to generation n 
there are Op{n) particles in all. Theorem 6 implies that at least a fraction 5 
of these are located at sites with (roughly) log re other particles. Thus, a 
substantial fraction of the particles fall in just Op{n/ \ogn) sites. This does 
not logically rule out the possibility that many more sites are occupied; 
however, it does suggest that the number f]„ of occupied sites is of order 
Op (re). This is consistent with the corresponding result for super-Brownian 
motion Xt, which states that for any t > 0, the random measure Xt is almost 
surely singular. Following is a sharp result about the number of occupied 
sites. 

Theorem 7. For two-dimensional nearest neighbor branching random 
walk, the number il„ of occupied sites is Op(re/logre) given the event Gn- 

Theorem 6 implies that the number of occupied sites must be of order at 
least re/ log re. Combining this with Theorem 7 we see that re/ log re is the 
true asymptotic rate. Revesz Revesz (1996) (Theorem 3 (ii)) asserts that a 
corresponding result is true for branching Brownian motion, but we believe 
that his proof has a serious gap. See section §7.2 for a detailed discussion. 

The next theorem partially quantifies the degree of particle clustering 
around a typical site. 

Theorem 8. Assume that d = 2. Let {in} any sequence of real num- 
bers such that lim.„£„ = oo and lim.„ log£„/log re = 0. Let Sn be the location 
of a typical particle, and let B{Sn',in) the ball of radius in centered at Sn- 
Then conditional on Gn, 

(A) the number of unoccupied sites in B{Sn',in) is op(£^), and 

(B) the number of particles in B{Sn',in) is of order Op(logre • i"^). 

Theorems 2 and 3 are proved in section §2, Theorem 4 in section §3, and 
Theorem 5 in section §4. Theorem 6 is proved in section §5, Theorem 8 in 
section §6, and Theorem 7 in section §7. For each of the last three theorems 
the calculations required for the proofs are considerably simpler in the special 
case of binary fission, where the offspring distribution Q is double- or-nothing 
- that is, Qo = Q2 = 1/2. In the interest of clarity, we shall give complete 
arguments only for this special case. These arguments (as should be evident) 
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can be extended to the general case of mean 1, finite variance offspring 
distributions. 

Fundamental to many of our arguments is the following elementary rela- 
tion between the expected number of particles at a site x in generation n 
and the n— step transition probabilities Pn{x) of the simple random walk: 

(1.9) EUn{x) = Pnix). 

This is easily proved by induction on n, by conditioning on the first gen- 
eration of the branching random walk. Here and throughout the paper, the 
term simple random walk is used for the symmetric nearest neighbor random 
walk on the lattice Z*^ with holding probability l/(2d+ 1) — that is, each in- 
crement is uniformly distributed on the set N" of 2d+l sites at distance < 1 
from the origin — and the notation Pn{x) is reserved for the probability 
that a simple random walk started at the origin finds its way to site x in n 
steps. We use the notation P" to denote the n— step transition probability 
kernel of simple random walk, that is, the nth iterate of the Markov operator 
P : £°°{Z'^) associated with the random walk. 

Notation. Following is a list of notation, in addition to that already estab- 
lished in equations (1.1), (1-2) and (1.9) above, that will be fixed throughout 
the paper: 

• M = {ej}^d<j<d is the set of sites at distance or 1 from the origin in 

• Q = {Qi}i>o is the offspring distribution, and Q* = {Q}} its ith con- 
volution power. 

• J-n is the cj— algebra generated by the random variables {Umix)}^^^d 

• A = 5/(47r) is the constant such that Pn{0) ~ A/n in dimension 2, see, 
e.g., P7.9 on Page 75 in Spitzer (1976) . 

In addition, we will follow the custom of writing / ~ to mean that the 
ratio f/g converges to 1, and f ^ g to mean that the ratio f/g remains 
bounded away from and oo. Throughout the paper, C,Ci,C' etc. denote 
generic constants whose values may change from line to line. Finally, we use 
a "local scoping rule" for notation: Any notation introduced in a proof is 
local to the proof, unless otherwise indicated. 

2. Proofs of Theorem 2 and 3. 

2.1. The case where the offspring distribution has finite moments. The 
proof of Theorem 2 will rely on the following estimates for the moments of 
the occupation statistics Unix). 
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Proposition 9. Suppose that the offspring distribution Q has finite ath 
moment for some integer a >2. 

(i) If d > 3, then 

sup V] ii;[/„(x)" < oo. 

n 

X 

(a) If d = 2, then there exist Ci,C2 < oo such that for all n, 

X 

Proof. We will use the following inequality: For alH > 2 and all bi > 0, 

(2-1) fe^^l <EEfE^^^-i{^>o>..>..>o}) ' 

\i=l / k=2 Vk \t=l / 

where Vk is the set of A;— tuples (ii, . . . ,ik) of distinct positive integers no 
greater than I. Inequality (2.1) is obviously true for I < a.To see that it holds 
for / > a, observe that, by the multinomial expansion, the left side of (2.1) is 
a sum of terms of the form t = (^.^ ^ j)^i^2 • • - ^i ^ where the exponents jj 
sum to a. Since at most a of these can be positive, and t vanishes if any of bi 
with exponent jj > is zero, the term t is included in the sum on the right 
side of (2.1). 

Next, by the Holder inequality, for each integer k >2 and all real numbers 
bi > 0, 



(fc \ ° k 
i=l / i=l 



This implies that if k independent branching random walks are started by 
particles ui, . . . ,Uk located at sites xi, . . . ,Xk respectively, and if U^^{x) is 
the number of the nth generation descendants at site x of the particle Ui, 
then 



)>0,...,U„'' {x)>0} 



.4 = 1 



(2.3) <A:--^ E ^^^n ix)r ■ n ^(^"' (^) > 0) 

i=l X j^i 

<A;" J^^C/„(x)" • (c— 



d 

n 
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Here we have used (2.2) in the first inequahty; the second inequahty follows 
by the local central limit theorem and the elementary observation that 

P{Un'{x) > 0) < EUn'{x) = Pn{x- X^). 

We are now prepared to estimate 'Y^^EUn{x)°^ . Conditioning on the first 
generation, we obtain 



yxT 



k 

Yl [Y^^nUix) ■ i{u:i,{x)>o,...,u:i,(x)>o} 



X 

a 

<Y,EUn-i{xr+Y.J2^ 

X k=2 X 

< ^ i^^„_r (xr . ( 1 + f; Q . f c-^ 

where Vk denotes the set of A;— tuples (ui, . . . ,Uk) of distinct particles in 
generation 1, and the first and second inequality hold by (2.1) and (2.3) 
respectively. Therefore, for all n, 

Eu^ixT < n ( 1 + E E • (^7=^) ) -E ^f^i 

Clearly, ^x^^^i-^)'^ — (2'^ + 1)EZ^ < oo. Furthermore, in dimensions 
d > 3, 

A.., /„ 1 



i=2 \ k=2 I ^ ^ \ 

\i=2 k=2 I ^ ^ \ 

=-p (^'EEq^I)^")' 

where C" < oo is independent of n; and in dimension d = 2, 



i=2 \ k=2 I 

< exp 



I ^ ^ i=2 k=3 I ^ ^ 1=2 ^ 



k-V 



<exp(C722" logn + Cs), 
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where C2 is a constant independent of both a and n, and C3 is a constant 
independent of n. □ 

Proof of Theorem 2 . By Kolmogorov's estimate (1.1), the probabil- 
ity that the process survives to time n is 0(l/n). By the Markov inequahty, 

P{Vn > Cn^^'^} < C^'^n~^EV^ < C^'^n^^ E^Un{xT , 

X 

and so the relation (1.3) follows from Proposition 9. □ 

Remark 10. Yaglom's limit theorem implies that, conditional on the 
event G.„, the number of particles at time n — 1 is Op{n). For each of 
these, there is a small chance that the number of offspring will exceed 
{2d + l)ni/(°+^), in which case K wih be at least ni/(°+=). If the tail of 
the offspring distribution decays like m~^°'~^^^ as m — t- 00, then the chance 
that one of the Op{n) particles in generation n — 1 will have more than 
{2d + l'jn^/('^+'^) offspring is of order one. Thus, the result in Theorem 2 is 
almost optimal. (This answers a question of Michael Stein.) 

2.2. The case where the offspring distribution has an exponentially de- 
caying tail. We begin with a stochastic comparison result for the random 
variables Un{x). First, observe that the law of the branching random walk 
(started by a single particle located at the origin) is invariant with respect 

to reflections in the coordinate axes, and so Un{x) = Un{x') for any two sites 
x' at corresponding positions of different orthants. Now define the usual 
partial order on the positive orthant U^: 

X ^ y if Xi < yi for all 1 < f < d. 

Lemma 11. If x ^y then Un{x) stochastically dominates Un{y); in par- 
ticular, Un{y) is stochastically dominated by Un{0) for every y ^U^. Conse- 
quently, if X ^ y, then for every n > 0, 

(2.4) Pn{x) > Pn{y) and 

(2.5) Un{x) > Un{y), 

where Un{x) := P{[7„(x) > 1} is the hitting probability function of the 
branching random walk. 

Remark 12. The relation (2.4), which follows from the stochastic dom- 
inance Un{x) >D Un{y) by taking expectations (recall the fundamental re- 
lation (1.9)), also follows more directly by the reflection principle for simple 
random walk. 
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Proof. Because the law of the branching random walk is invariant with 
respect to permutations of the coordinates, we may assume, without loss of 
generality, that y = x + ei, where ei = (1, 0, . . . , 0). Denote by L and L' the 
hyperplanes 



observe that y is the reflection of x in L' . We shall deflne a particle system 
with particles of three colors — red, blue, and green — in such a way that 

(a) the subpopulation of all red and blue particles follows the law of the 
branching random walk started by one (red) particle at the origin; 

(b) the subpopulation of all red and green particles follows the same law; 

(c) there are no red particles to the right of the hyperplane L'; and 

(d) at each time, the green and blue particles are paired (bijectively) in 
such a way that the green and blue particles in any pair are at sym- 
metric locations on opposite sides of the hyperplane L' . 

This will prove that Un{x) >t> Un{y) for each n, by the following reason- 
ing: First, the distribution of Un{x) coincides with the distribution of the 
total number of red and blue particles at location x and time n, by (a). Sec- 
ond, the number of blue particles at x equals the number of green particles 
at y, by (d), since x and y are at symmetric locations on opposite sides of 
the hyperplane V . Third, the number of green particles at y has the same 
distribution as Un{y), by (b) and (c). 

The particle system is constructed as follows. To start, color the initial 
particle at the origin red. Offspring of blue and green particles will always 
have the same color as their parents, and each blue particle b will always 
be paired with a green particle g located at the mirror image (relative to 
reflection in the hyperplane L') of the site of b. Offspring of red particles will 
be red except possibly when the parent red particle is located at a site on 
the hyperplane L. In this case — say, for definiteness, that the red parent 
particle ^ is at site z £ L — each offspring particle ^ first makes a jump 
according to the law of the nearest neighbor random walk, and then chooses 
a color as follows: (a) If the jump is to a site z' ^ zto the left of hyperplane L' 
then C becomes red; and (b) If the jump is either to the same site z as the 
parent or to its mirror image z* on the right of L' then ^ chooses randomly 
between blue and green. In case (b) the offspring particle generates a 
doppelganger (mirror particle) ^' of the opposite color at the reflected site 
on the other side of L' . Note the distribution of the position of is the same 
as that of The particle C generates an offspring branching random walk 



L = {z€W 
L' = {z£ M' 



d 



z\ = xi} and 
zi = xi + 1/2}; 



d 
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with all particles having the same color as the mirror image Q(^i of 
relative to L' (with particles colored oppositely) is attached to Note that 
Q^f is itself a branching random walk started at the location of d^' , by the 
symmetry of the nearest neighbor random walk. 

Properties (a)-(d) above are now readily apparent. Property (c) holds 
because, by construction, children of red particles on L that jump across L' 
are either green or blue, and offspring of blue and green particles are either 
blue or green. Property (d) is inherent in the construction. Finally, (a) and 
(b) follow from the blue/green symmetry of the reproduction law for red 
particles located at sites on L. □ 

Proposition 13. Assume that the offspring distribution Q has finite 
moment generating function in some neighborhood of the origin. Then in 
dimensions d> 3, there exist 5d > and C > such that for any 9 € [0, 6ii], 
all X and all n > 1, 

(2.6) Ee^p{eUn{x)} - 1 < CPn{x)e. 

In dimension d = 2, there exist ^2 > and C > such that for any E [0, ^2], 
all x G and all n>l, 

(2.7) Eexp{^C/„(x)/logn} - 1 < CP„(x)^/ log n. 

Proof. Let ^{z) = Xli^i Qi^^ probability generating function 

of Q. By hypothesis, $(2) is finite and analytic in a neighborhood of the 
closed disc \z\ < for some 5 > 0, and since the variance of Q is strictly 
positive, ^{z) is strictly convex on [0, e^]. Moreover, $'(1) = 1, because the 
offspring distribution has mean 1. 

Define 

Gn{x) = Gn{x;e) = Eexp{eUn{x)) - 1. 

Clearly, Gn{x;6) — )• as ^ — )• 0. Moreover, by Lemma 11, for each value of 
6 > the function Gn{x) is maximal at x = 0. Since the random variables 
Ui{x) are zero except for x G AA, and have the same distribution for x € N , 
the function Gi{x) is, for any fixed 6, a scalar multiple of the uniform dis- 
tribution Pi on J\f . Conditioning on the first generation of the branching 
random walk shows that 

(2.8) G„+i(x) + l = $(PG„(x) + l) 

where P is the one-step Markov operator for the simple random walk, that is, 
P/(x) = Ef{x + Y) where Y is uniformly distributed on N . Since <^'(1) = 1 
and $(^) is strictly convex for z G [0, e^], equation (2.8) implies that 

(2.9) Gn+i{x) < PG„(x)$'(l + PG„(x)). 
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Unfortunately, both relations (2.8) and (2.9) are nonlinear in G„. For this 
reason, we introduce dominating functions Hn{x) = Hn{x;6) that satisfy 
corresponding linear relations: Set Hi[x) = G'i(x), and define Hn inductively 

by 

(2.10) Hn+i{x) = rHn{x)<t>'{l + F„(0)). 

Note that Hn+i may take the value +oo if Hn{0) exceeds the radius of 
convergence of ^. Since Gn{x) < G„(0), the inequality (2.9) implies that 
H2 > G2, and so by induction that Hn > Gn for all n > 1. Thus, to prove 
inequalities (2.6) and (2.7) it suffices to prove analogous inequalities for the 
functions Hn{x;9). 

The advantage of working with the functions Hn is that the linear relation 

(2.10) can be iterated. In general, if functions / and g satisfy g = aPf for 
some scalar a, then Fg = aP^/. Employing this identity in equation (2.10) 
and iterating yields 

n-l 

Hn{x) = ¥^-'Hiix) H <!>'(! + Hj{0)). 

i=i 

Because the function Hi = Gi is itself a scalar multiple of Pi , it follows that 

(2.11) Hn{x; 9) = Pn{x)Hi{0; e){2d + 1) H ^'(^ + ^))- 

i=i 

Since ^"'(1) = 1, the factors in the product are well-approximated by (1 + 
^"{l)Hj{0;9)) as long as Hj{0;6) remains small. In particular, for suitable 
constants C < 00 and e > 0, if Hj{0; 6) < e for all j < n — 1 then 

(2.12) Hn{x;e) < {2d+l)Pn{x)Hi{0;e)l[{l + CH,{0;9)), 
equivalently. 

The large-n behavior of the products on the right side of (2.12) will depend 
on whether or not the sequence Pn{0) is summable, that is, on whether or 
not the simple random walk is transient. There are two cases to consider: 

Dimensions d > 3: In dimensions d > 3, the return probabilities Pn(0) are 
summable. Moreover, when 9 > is small, the factor {2d + l)Hi{0; 0) on the 
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right side of (2.13) is also small, because Hi = Gi is a continuous function 
of 9 that takes the value at = 0. Hence, by choosing 9 small we can 
make the sum over n of the quantities on the right side of inequality (2.13) 
arbitrarily small. Now the fraction on the left side of (2.13) is the nth term 
of the telescoping series 

(2.14) C~'Y1 



^nU(^ + CH,{0;9)) UU(^ + CH,iO;9))J' 
consequently, (2.13) implies that for all sufficiently small 9 > the products 

n 

ll{l + CH,{0;9)) 
i=i 

remain bounded for large n, and for small 9 remain close to 0. It now follows 
by (2.12) that for a suitable constant C' < oo and all small 9 the functions 
Hn{x;9) are all finite, and satisfy 

Hnix;9) < C'Pn{x)Hi{0;9). 

Finally, the differentiability of Hi{0;9) in 9 guarantees that Hi{0;9) < C9 
for an appropriate constant C < oo for all small 9. This proves (2.6). 

Dimension d = 2: It is still the case that the fraction on the left side of 
(2.13) is the nth term of the telescoping series (2.14), but since X]-Pn(0) 
diverges, this no longer implies that the products on the right side of (2.12) 
remain bounded. However, the local central limit theorem gives an explicit 
estimate for the partial sums of the return probabilities: in particular, for 
some C' >A = 5/(47r), 

n 

^ Pj (0) < C log n for all n > 2. 

Consequently, substituting 9/\ogn for 9 in inequality (2.13) and summing 
gives 

n 

1 - J\{l + CHj{Q;9/\ogn))-^ < C"9. 
i=i 

This in turn implies that 

n 

+ CHjiO; 9/ log n)) < 1/(1 - C"9). 

i=i 

Using this upper bound for the product on the right side of (2.12) and using 
the bound Hi{0; 9/logn) < C9/logn for smah 9 yields (2.7). □ 



imsart-aop ver. 2009/08/13 file: BRWpublish. tex date: April 8, 2010 



14 



STEVEN P. LALLEY AND XINGHUA ZHENG 



Remark 14. In dimensions d > 3, the conclusion (2.6) cannot be ex- 
tended to all ^ > 0, even for the double-or-nothing case. In fact, for suffi- 
ciently large 9, the sums YlxeZ'' (-^ 6xp{0C/„(x)} — 1) are not bounded in n. 

Remark 15. In dimension d = 2, the relation (2.7) does not hold for 
large 9. See Remark 26 below. 

We are now prepared to prove Theorem 3. In fact, we will establish the 
following stronger result: 

Corollary 16. Under the hypotheses of Proposition 13, with the same 
notations. 



(i) If d> 3, then for all 9 < 6^, 

log n 



o 



1 



n 



In particular, conditional on Gn, Vn = Op{logn). 

(ii) If d = 2, then for all 9 < 62, 

(logn)^ 



P Vr,> 



Gn 



o 



1 



nS2/e - 1 ; ■ 

In particular, conditional on Gn, Vn = Op((Iogn)^). 



Proof. We will prove this only for dimensions d >3; the dimension d = 2 
case can be handled similarly. By Markov's inequality, 



P\Vn> 



logn 



Gn < 



1 



exp{6d/9 ■ logn] 



{exp{6dUn{x)) ■ l{U„(x)>0}\ Gn) 



^ i^ nsje_i • (exp((5d[/„(x)) • l{c/„(^)>o})^ • 



For any random variable X > 0, 



Eexp{X) = Sexp(X) • l|;^>o} + ^exp(X) • l^x=o} 
= EeMX)-l{x>o} + P{X = 0) 
= EeMX) ■ l{x>o} + l-P{X>0). 
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{exp{6dUnix)) ■ l{u„{x)>o}) 

X 

= Y,{Eexp{6dUn{x)) -l) + Yl P{Un{x) > 0) 

X X 

<^iEexp{6dUnix))-l) + l, 

X 

SO by Proposition 13, 

E {exp{6dUn{x)) • l{(7,^(^)>o}) < C for all n > 1. 

X 

The conclusion follows. □ 

3. Proof of Theorem 4. The proof uses the following elementary lemma, 
whose proof is left to the reader. 

Lemma 17. Suppose that on some probability space (Q, T ^ P) there are 
two events Ei,E2 such that 

(31) PjE.AE,) ^ 

where E1AE2 is the symmetric difference of Ei and E2- Then 

(3.2) \\P{-\Ei) - P{-\E2)\\tv < 2e, 

where P{-\Ei) denotes the conditional probability measure given the event Ei 
and II • \ \tv denotes the total variation distance. 

Lemma 17 will allow us to replace the event of conditioning G„ in Theo- 
rems 4 and 5 by asymptotically equivalent events of the form 

(3.3) Hn = {Z^(n) > nsn}. 

Lemma 18. Let m[n) < n be integers and En > real numbers such that 
m{n)/n — t- 1 and e„ — t- as n —t- 00. Then 

(3.4) lim — , , = 0. 

n^oo F(Gn) 
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Proof. This is an easy consequence of Kolmogorov's estimate (1.1) and 
Yaglom's theorem for critical Gallon- Watson processes. Let Kn = > 1}. 

Clearly, Hn C K^, and so P{Kn \ H^) = 1- On the other hand, Yaglom's 
theorem implies that P{Hn \ Kn) — )• 1, since m[n)/n — )■ 1. Consequently, 

(3.5) lim — — — — = 0. 

n-s>oo J-'[Kn) 

A similar argument shows that the symmetric difference KnAGn is an 
asymptotically negligible part of J^„. Obviously, G„ C Kn, so P{Kn \ Gn) = 1- 
Yaglom's theorem implies that for any J > there exists a > such that 

P{Zm{n) > an I Kn) >l-6. 

But on the event {Zjn[n) > Q;n} the event G„ of survival to generation n 
is nearly certain for large n, because the ^^{n) particles in generation m{n) 
initiate independent Galton- Watson processes, each of which survives to gen- 
eration n with probability ~ 2/ {n—m{n))a'^ , by Kolmogorov's estimate (1.1). 
Hence, 

P{Gn \Kn)>l-26 

for large n. Since 5 > is arbitrary, it follows that P{Gn \ Kn) — )• 1- By 
Lemma 17 and (3.5) we get 

(3.6) P{Gn\Hn)^l. 

Furthermore, since Gn C Kn-, 

n-5>oo P[Kn) 

By Lemma 17, this implies that conditioning on Gn is asymptotically equiva- 
lent to conditioning on Kn, and so the difference P{Hn \ Kn) — P{Hn \ Gn) — >■ 0. 
But we have seen that P{Hn \ Kn) — )• 1, hence P{Hn \ Gn) — ?• 1- This, along 
with (3.6), implies (3.4). □ 

Proof of Theorem 4. The offspring distribution is non-degenerate, so 
there exists Iq > 1 such that Qi^ > 0. Let p = Qi^^ ■ {l/{2d+ 1))'" be the 
probability that the initial particle produces Iq offspring and these offspring 
all stay at the origin. Then for all A; G N, 

-P(C/fc(0) > Iq) > p ■ p'" • p'" . . .p'o~' > p^o/ilo-^)_ 
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Our objective is to show that for some 6 > 0, P{Vn > 51ogn|Z„ > 0) — t- 1. 
By Lemmas 17 and 18, this will follow if we can show that for some m(n) < n 
with m{n)/n — t- 1 and some e„ — t- 0, the probability 

P{Vn > 5 log > nsn) 1. 

To do so, for (5 > to be determined later, and all n big enough, define k 
such that lodlogn > Iq > 6logn, and m(n) =n — k. Then m{n)/n — )• 1. Fix 
a sequence e„ = 0(1/ log n); then 



P K > elogn 



>i-(i-m(o)>zS) 

>l-(l-/V('o-l))'"" 



n 



>1 -exp (e„n (-p'o/('o-i) 



(3.7) 

>1 - exp (^-En 



^p'o51ogn/(io-l) 



=1 - exp (^-e„ni+'«'^'°s?'/(^o-i) 
provided that 6 < {Iq — l)/( — logp). □ 
4. Proof of Theorem 5 . 



4.1. Strategy. By Lemmas 17 and 18, the difference between conditioning 
on the event Gn = {Zn > 0} and conditioning on the event Hn := {Z^^(^n) ^ 
nSn} is asymptotically negligible if ?n(n)/n — )■ 1 and e„ — )■ 0. Thus, it suffices 
to prove the weak convergence of the conditional distributions in (1.7) when 
the conditioning event is Hn rather than The advantage of this is that, 
conditional on the state of the branching random walk at time m(n), the next 
n — m{n) generations are gotten by running independent branching random 
walks for time n — m(n) starting from the locations of the particles in gener- 
ation m(n). The argument will hinge on showing that if m(n) < n is chosen 
appropriately then these independent branching random walks will not over- 
lap much at time n, and so the total number M„(j) of multiplicity-j sites 
will be, approximately, the sum of .^m.(n) independent copies of M„„^(„)(j). 

4.2. Overlapping. 

Lemma 19. Suppose that a critical branching random walk starts at time 
with two particles u,v located at sites Xu,Xy S Z,'^, respectively. Let Dn{u,v) 
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be the number of particles in generation n located at sites with descendants 
of both u and v. Then there exists C > such that for all generations n > 1, 

(4.1) EDn{u,v) < 2P2n{x^-Xu) < C (l/V^)^ 

Proof. Denote by Un{x) the number of descendants of particle ^ at site x 
in generation n. Since the progeny of particles u and v make up mutually 
independent branching random walks, the random variables U^{x) and U^{x) 
are independent. But 

EDn{u, v)=EY^ (Unix) + U^{x)) l{U}t{x)>l} ^{U-{x)>l} 

= 2 ^ EU^{x) l{c/^(x)>i} 

^ 2 ^ ^ Pn(x X^^P^iix X-ij) 

x&,<i 
— "^PiniXy Xu) 

<c{i/v^'. 

□ 

Corollary 20. Let Yn-m be the number of particles in generation n 
located at sites with descendants of at least two distinct particles of generation 
m < n. Then 

(4.2) E{Yn-m I J'm) < CZi/{n - mfl'^. 
4.3. Convergence of means. 

Proposition 21. In dimensions d>3, 

(4.3) lim£^M„(j) = Hj exists for every j > 1, and 

n 

oo 

(4.4) Ej'-^^ = 1- 
i=i 

Proof. The random variable Mn{j) counts the number of multiplicity— j 
sites in generation n. The particles at such a site will either all be descendants 
of a common first-generation particle or not; hence, by conditioning on the 
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first generation of tlie branching random walk we may decompose Mn+i{j) 
as follows: 

Zi 

(4.5) M„+i(j) = Miij) + A„+i(j) - S„+i(i) 

1=1 

where (a) the random variables {M^(j)}j<^^ are independent copies of Mn{j); 
(b) the error term An+i{j) is the number of multiplicity— j sites at time n+1 
with descendants of different particles in generation 1; and (c) the correction 
Bn+i{j) equals 

# particles in generation 1 

xeMn+i(j+) 

with exactly j descendants at x in generation (n + 1), 

where A^n+i(j+) is the set of sites with (j + l) or more particles in generation 
(n + 1). Obviously, A„_|_i(l) = 0, because a site with only one particle cannot 
have descendants of distinct first generation particles, and so it follows that 
EMn+i{l) < EMn{l). This implies that lim„SM„(l) exists. 

To see that Yvoin^o^ EMnij) exists for j > 2, observe that both An+i{j) 
and Bn+i{j) are bounded by the number of (n + 1)— th generation particles 
at sites with descendants of different particles of generation 1. Hence, by 
Lemma 19, writing Z{1) = Z\ for the first generation of the branching 
process, 

(4.6) ^(A„+i(i) + B„+i(i)) <2S Dn{u,v) 

1=2 ^ ^ 

< C'n-''/\ 

for some C" < oo, because the offspring distribution has finite second mo- 
ment. Consequently, by equation (4.5), 

\EMn+l{j) - EMn{j)\ = 0(71-''/^). 

Since the sequence vT'^l'^ is summable for (i > 3, the sequence {£'M„(j)}„>i 
must converge. This proves the convergence of means (4.3). 

Clearly, for each n > 1 it is the case that jEMn{j) = EZn = 1. 
Hence, to prove the equation (4.4), it suffices to show that for every e > 
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there exists an integer k = k[e) such that for all n > 1, 

oo 

(4.7) EYn{k)<e where y„(A;) = j • M„(i) 

j=k 

is the number of particles in generation n located at sites with at least 
{k — 1) other particles. Since Yn{k) < > /c}, and since EZ^ = 1, it is 

certainly the case that for any fixed n > 1 and e > there exists k = k{n; e) 
so that inequality (4.7) holds; the problem is to prove that k{e) can be chosen 
independently of n. By the same reasoning as in relation (4.5) above, for all 
n,k > 1, 

(4.8) Yn+i{k)= Y:ik) + Cn+i{k) 

where the random variables Yj^{k) are independent copies of Yn{k) and the 
error term Cn+iik) is bounded by the total number of particles in generation 
n+1 at sites with descendants of at least two distinct particles in Z{1). Since 
EZi = 1, the decomposition (4.8) implies that 

\EYn+i{k) - EYn{k)\ < ECn+i{k). 

But by the same logic as in relation (4.6) above, there exists C < oo inde- 
pendent of k and n such that ECn+i{k) < C'n~'^/'^ for all 77,, > 1. It follows 
that for sufficiently large ?i(e) and all /c > 1, 

oo 

ECn+i{k)<e. 

n=n(e) 

Thus, if for some k > \ and n = n{e) the inequality (4.7) holds, then 
EYn{k) < 2e for all n > n{e). This proves (4.4). □ 

Remark 22. Since the error term Cn+i{k) in equation (4.8) is nonnega- 
tive, the expectations EYn{k) are nondecreasing in 77. Because the offspring 
distribution is nondegenerate, for every A; > 1 there exists n > 1 such that 
Yn{k) > 1 with positive probability, which forces EYn{k) > 0. Therefore, 
there are infinitely many integers j > 1 such that Kj > 0. 

4.4. Conditional weak convergence: Proof of Theorem 5 . In view of Kol- 
mogorov's estimate (1.1), the inequality (4.7) can be rewritten as 

E \YjMn{j) gA < Cne 
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for some constant C < oo not depending on n. Since r2„, = Mn{j), it 
follows that to prove Theorem 5 it suffices to prove that for any finite A; > 1, 



(4.9) 




j<k> 



where Y is exponentially distributed with mean For this, we will use 

Yaglom's theorem, the convergence of moments (4.3), and a crude bound on 
the variance of Mn{j)- 

(4.10) Var(M„(j)) < EZ^ = 1 + na^. 

Fix 1 < m < n, and for each particle u € let M^_^(j) be the number 
of sites that have exactly j descendants of particle u in generation n. The 
random variables M^_^{j) are, conditional on J^rn, independent copies of 
Mn-m{j)- Now Mn{j) decomposes as 

(4.11) M„(j)= Yl M„"_„,(j) + i?„;™,^M;^^(i) + i?„;^ 

where the remainder Rn;m is bounded, in absolute value, by the number of 
particles in generation n located at sites with descendants of at least two 
distinct particles of generation m < n. By Corollary 20, 

(4.12) E{\Rn-m.\ I -F™) < CZl/{n - mfl''. 

By Yaglom's theorem, the conditional distribution of Zm/m given the event 
Gm of survival to generation m converges to the exponential distribution 
with mean 2/o"^; thus, if m = m{n) is chosen so that m/n — )• 1 and n — m > 
^2/{d-e) £qj, gQj^g e > 0, then the bound in (4.12) will be of order op{n). In 
view of (4.11) and Lemmas 17 and 18, it follows that to prove (4.9) it suffices 
to prove the corresponding statement in which the random variables M„(j) 
are replaced by the approximations M*.^(j) in (4.11), and the conditioning 
events Gn are replaced by the events Hn = {Zm > ^nn}- But this follows 
routinely by first and second moment estimates: if the scalars e„ are chosen 
so that En — ?• but ne„/(n — m) — )• oo, then by relation (4.3) and the variance 
bound (4.10), 



u£Z{m) 
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and 

Var I Yl ^n-mU) J'^n | < + (n - m)a^)/Z™ 0. 

\ u£Z{m) 

Chebychev's inequality now implies that the conditional distribution of 
M*.^{j)/Zm given Hn is concentrated in a vanishingly small neighborhood 
of Kj as n — )• oo. Since the conditional distribution of Zm/n given Hri con- 
verges to the exponential distribution with mean by Yaglom's Theorem 
and Lemmas 17 and 18, the desired result follows. □ 

5. Typical Sites in Dimension 2: Proof of Theorem 6. 

5.1. Embedded Galton-Watson tree. For simplicity we shall consider only 
the binary case, that is, the special case where the offspring distribution is 
the double-or-nothing distribution Qo = Q2 = 1/2- The arguments can all be 
easily adapted to the general case, at the expense of notational complexity. 

We begin with the simple observation that the branching random walk 
can be constructed by first generating a Galton-Watson tree r according to 
the given offspring distribution, then independently attaching to the edges 
of this tree random steps, distributed uniformly on the set N of nearest 
neighbors of the origin. The vertices of r at height n represent the particles 
of generation n; the location in 1? of a particle a of the nth generation is 
obtained by summing the random steps on the edges of the path in r leading 
from the root to a. Henceforth we will distinguish between the underlying 
Galton-Watson tree r and the marked tree r* obtained by attaching step 
variables to the edges of r. Observe that the conditional distribution of the 
marks of r* given the tree r is the product uniform measure on Af^^'^^ , where 
£{t) denotes the set of edges of r. 

A typical particle of the nth generation in a branching random walk con- 
ditioned to survive to the nth generation can be obtained by first choos- 
ing a tree r randomly according to the conditional distribution Fn of the 
Galton-Watson tree given the event of survival to generation n, then ran- 
domly selecting one of the Z„, > 1 vertices at height n. For this random 
choice we assume that the underlying probability space supports a uniform- 
[0, 1] random variable 7 independent of all other random variables used in 
the construction of the branching random walk. Since this procedure does 
not use information about the step variables attached to the edges of the 
tree, it follows directly that the trajectory of the typical particle, conditional 
on the underlying Galton-Watson tree, is a simple random walk started at 
the origin. 
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5.2. Reduction to the Size-Biased Case. The strategy of the proof of The- 
orem 6 win be based on a change of measure. Denote by Ph = Pj^ the 
probabihty measure that is absolutely continuous relative to P with Radon- 
Nikodym derivative 

(5.1) = Z^. 

The measure P^ so defined is a probability measure, because EZn = 1. Call 
it the size-biased measure. In the arguments below the value of n will be 
fixed, so we will generally omit the dependence of the measure on n and 
write Ph = Pj^- Because the Radon-Nikodym derivative depends only on 
the underlying Galton- Watson tree r, which under P is independent of the 
marks, it follows that the conditional distribution under Ph of the marks 
given the tree r is the same as under P. Thus, to construct a version of 
the marked tree r* under Ph, one may first build a size-biased version of 
the underlying Galton- Watson tree, then attach edge marks independently 
according to the (product) uniform distribution on M. Henceforth we will 
call such a marked tree a size-biased marked tree or a size-biased branching 
random walk. 

Observe that Ph is also absolutely continuous relative to the conditional 
distribution P* of P given the event Gn of survival to generation n; the 
Radon-Nikodym derivative is 

(5.2) ^ = ZnTTn 

where vr^ = P{Gn) ~ (2/ncr^). By Yaglom's theorem, under P* = P{- \ Gn) 
the distribution of dP^ / dP* converges in law to the unit exponential distri- 
bution. This implies the following. 

Lemma 23. To prove Theorem 6 it suffices to prove the analogous state- 
ments for the measure Ph, that is, to prove that (i) for each e > there 
exists K < oo such that 

(5.3) PH{Tn> K log n} <e- 

and (a) for all sufficiently small e > there exists 5 > such that for all 
large n, 

(5.4) PH{Tn>elogn}>5. 

Proof. This is a direct consequence of the fact that the Radon-Nikodym 
derivatives dPn/dP* converge in law under P* as n — )■ oo, because this 
implies that the Radon-Nikodym derivatives dP* /dPn converge in law un- 
der Ph- □ 
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5.3. Structure of the size-biased process. The size-biased measure Ph 
on marked trees is especially well-suited to studying typical points, and 
has been used by a number of authors (see Lyons, Pemantle, and Peres 
(1995) and the references therein) for similar purposes. Consider first the 
distribution of the unmarked genealogical tree r under Ph- According to 
Lyons, Pemantle, and Peres (1995), a version of this random tree can be 
obtained by running a certain Galton-Watson process with immigration. In 
the case of the double-or-nothing offspring distribution, the nature of this 
process is especially simple: 

Recipe SB: Each generation j has a single distinguished particle Vj which 
gives rise to two particles in generation j ' + 1, one the distinguished parti- 
cle Uj+i, the other an undistinguished particle. All undistinguished particles 
reproduce according to the double-or-nothing law. For each n, the distin- 
guished particle v„ is uniformly distributed on the particles in generation n. 

□ 

Thus, a version of the size-biased branching random walk, together with 
a randomly chosen point f.„ of the nth generation, can be built by attaching 
independent step random variables to the edges of the random tree built 
according to Recipe SB. Equivalently, this process can be constructed using 
three independent sequences of auxiliary random variables: 

(Tq) {«S'„}„>o is a simple random walk in 1? with initial point Sq = 0; 

(^f*) {^i}«>0 ^-re independent and uniformly distributed on A/"; 

(Tc) {[/^(a;)}j>o are independent copies of the branching random walk {?7„(a;)} 

run according to the law P; and 
(Td) Bq ~ Bernoulh(l/(2d + 1)). 

(We emphasize that the auxiliary branching random walks {C/^(a;)}j>i are 
run according to the original probability measure P, not the size-biased 
measure Ph-) The size-biased branching random walk is obtained by letting 
the "typical" particle follow the trajectory Sj, then attaching an additional 
particle to each point {j,Sj) visited by the typical particle, letting it make 
a step to Sj + S^j, and then attaching the jth copy of the branching random 
walk to this particle. 

Corollary 24. The distribution of Tn under the size-biased measure 
Ph is the same as the distribution under P of the random variable 

n-2 

(5.5) T* = i + B^ + Y, U^-j-iiSn - S, - e,). 

j=0 
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□ 

The Bernoulli random variable Bq accounts for the possibility that the 
sibling of the typical particle jumps to the same site as the typical particle. 

Reversing the random walk will not affect the distribution of the random 
variable Tn, since the random walk is independent of all other component 
variables of the representation (5.5), nor will reversing the indices of the 
auxiliary branching random walks . Thus, the following random variable 
has the same distribution as that given by (5.5): 

n 

(5.6) rr = 1 + i?o + E + O-i)- 

i=2 

5.4. Variances of the occupation random variables. Next we focus on the 
distribution of the random variable T** defined by (5.6). To obtain con- 
centration results for this distribution, we will need bounds on the second 
moments of the random variables [/„(x); for this, we use an exact formula 
for the second moment of f/„(x), valid in all dimensions: 

Proposition 25. 

n-l 

(5.7) EUnixf = Pn{x) +a^Y.Yl P^(^)Pn-iix " z), 

i=0 z 

Proof. This is a special case of equation (81) in Lalley (2009), which 
gives the mth moment for all integers m > 1. In the case m = 2, a simple 
proof can be given by conditioning on the first generation of the branching 
random walk. Set fn{x) = EUn{x)'^ and gn{x) = Pn{x)'^; then conditioning 
on generation 1 gives 

fn{x) = ¥U-i{x)+a^gr,{x). 

Since the operator P is linear, this relation can be iterated n — 1 times, 
yielding 

fn{x) = P"-Vl(^) + a^Y.^'gn^i{x). 

i=0 

This is equivalent to the identity (5.7). □ 

Remark 26. If the offspring distribution has an exponentially decaying 
tail, then one can deduce from (2.7) that EUn{x)'^ < Clogn/9. However, 
formula (5.7) implies that X^^, -E'C^n(a;)^ grows at rate logn, so (2.7) cannot 
hold for large 6. 
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5.5. Mean and variance estimates for T** . The sum in the representation 
(5.6) can be decomposed as r„ + A^, where 

n 

(5.8) Tn: = Y.P^{Si) and 

1=2 
n 

(5.9) A„: = ^X,_i with Xi^i:=Uizl{S, + ^i^i)-P,{Si). 

1=2 

Lemma 27. Let Sn be simple random walk in 7?, and let r„ be defined 
by (5.8). Then 

ETn A 

(5.10) lim = — and 

n-s>oologn 2 

(5.11) hm Var f ) = 0. 

n^oo \^Og n J 

Recah that A = 5/(47r) is the constant such that Pn{0) ~ A/n. 

Proof. By the symmetry of the simple random walk, EPi{Si) = P2i{^) ~ 
A/{2i), and so the first convergence (5.10) follows routinely. To estimate the 
variance, first observe that 

n 

EVl = 2^EP,{S,)P,{S,) + ^EP,{S,f 

i<j i=2 

(5.12) =2Y,EP^iS^)PJ{SJ) + 0{l). 

i<j 

The second equation follows from the local central limit theorem in d = 2, 
which guarantees that Pi{z) < C/i for some constant C < oo independent of 
i and z. Next, observe that for i < j, by the symmetry of the random walk 
and the fact that Pi{z) is maximal at z = (Lemma 11) 

(5.13) EPi{Si)Pj{Sj) = E{E{P,{Si)PjiSj)\S,)) 

= EPiiSi) Pj{S^+x)Pj^^{x) 

= EPi{Si)P2j-i{Si) 

< EPi{Si)P2j-^{0) 
= P2^(0)P2i-^(0). 
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Substituting this bound in (5.12) and applying the local central limit theorem 
(in the form Pn{0) ~ A/n) yields 

n 

Y,m{s.)p,{Sj) < 5] J]P2^(o)P2,-^(o) 

i<j j=2 i<j 

n 

< 2 ^ ^ A'^/{2i{2j - i)) + error 
j=2 i<j 

^\ 2 
~ — log n + error, 
4 

where the error is of smaller order of magnitude. Together with (5.12) and 
(5.10), this shows that 

Var (r„) = EVl - {EFnf = o{lognf. 

□ 

Lemma 28. Let Sn, U^{x), and be independent sequences of random 
variables satisfying the hypotheses (T^) — (T^) of section §5.3. If An and Xi 
are defined as in equation (5.9), then 

(5.14) EXi = and EXiXj = for all i / j. 

Consequently, 



(5.15) EAn = and lim Var 



n— >oo 



log n 



Proof. To show that EAn = it suffices to show that EXi = 0. This 
follows from the fundamental relation (1.9) by conditioning on Sj+i and ^j: 

EXi = EE{Ui{Si+i + ii) I Si+i, ii) - EPi+,{Si+i) 
= EPi{Si+i + ii) - EPi+i{S,+i) = 0. 

Now consider the covariances EXiXj. To compute these expectations for 
i < j, condition on the random variables /Sj+i, -Sj+i, {Ui (x)} , and (but 
not ^j), and use the fundamental identity (1.9): This implies that EUj{x + 
^j) = Pj_|_i(x) for each x G Z^, and so 

EXiXj = EEi^XiXj I •) 

= EXiE{U^^{Sj+i + i,) - Pj+i{S,+i) I •) 
= EXi -0 = 0. 
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It follows that the variance of the sum A„ is the sum of the variances of the 
increments Xi, and so 

n n 

Var(A„) = Y,EXl, = Y,EXU 

i=2 i=2 

n 

i=2 

= Y^EUizl{S. + i^-if + 0{l). 

Now by the second moment formula (5.7), 

EUizliSi + ii.if 

=E {pi-i{S, + ^i-i) Pjiz?Pi-j-i{Si + ^i-i - z) 

i-l 

j = l z 

=P2.(0) + 5]^P,(z)2.P2.-,(^). 

j=l z 

The first term is of order 0{l/i). To estimate the second, observe that by 
the local central limit theorem, for large j, 

where [•] denotes integer part and the relation holds uniformly for 1^1 < C-y/J. 
Consequently, for large i, 

j=l z j=l z 

A A 



E 



log i 
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Summing from i = 1 to n a shows that Var(A„) ~ (^^/8)log^n. This 
proves (5.15). □ 

5.6. Proof of Theorem 6: Binary fission case. By Lemma 23, it suffices 
to prove assertions (5.3)-(5.4). By Corollary 24, the distribution of T„ under 
the size-biased measure Ph is identical to the distribution of the random 
variable T** := 1 + under P, where T** is defined by (5.6). Finally, 

by Lemmas 27 and 28 (note that i^A^r^ = 0), 

A ~ 

ETn ~ — log n and Var(T„) — log^ n. 

Z o 

The first of these implies, by the Markov inequality, that T** = Op(logn). 
This proves the first assertion (i) of Lemma 23. The second assertion (ii) is a 
consequence of the following elementary lemma (see, e.g., Lawler and Limic 
(2007), Lemma 12.6.1). 

Lemma 29. If X is a nonnegative random variable with positive, finite 
second moment, then for any a E [0, 1], 

(5.16) P{X >aEX} > {I- afiEXf/EX"^. 

□ 

6. Clustering in Dimension 2: Proof of Theorem 8. 

6.1. Occupied sites in the ball B(Sn',^n)- We consider only the case of 
binary fission. The proof of Theorem 8 in this case, like that of Theorem 6, 
is based on the change of measure strategy outlined in section 5.2. In partic- 
ular, we shall prove the corresponding assertions to statements (A)- (B) of 
Theorem 8 for the size-biased process of section 5.3. Thus, assume through- 
out this section that the random variables Sj, Uj'^, and are as in (T^), 
(Tfe), (Tc) of section 5.3. Recall that the size-biased branching random walk 
is obtained by letting the "typical" particle follow the trajectory Sj, then 
attaching an additional particle to each point {j,Sj) visited by the typical 
particle, letting it make a step to Sj + ^j, and then attaching the jth copy 
of the branching random walk to this particle. To prove Theorem 8 it 
suffices to prove the following proposition. 

Proposition 30. Let {in} be any sequence of real numbers such that 
lim„^„ = oo and limjilog^n/logn = 0. Let B{Sn',in) the ball of radius in 
centered at Sn- Then for the size-biased branching random walk, 
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(A) the number of unoccupied sites in B{Sn',^n) is op{£'^), and 

(B) the number of particles in B{Sn'An) is of order Opilogn • 

The construction of section 5.3 shows (cf. formulas (5.5) and (5.6)) that 
the number of particles at location Sn + x in the nth generation of the size- 
biased branching random walk is distributed as 

n.-l 

(6.1) U:*{Sn + x) := 6o{x) + Bo ■ l{|,|<i} + ^ U^{Sj+i +x + ^,). 

i=i 

6.2. Vacant Sites: Proof of Theorem 8 (A). The representation (6.1) im- 
plies that the probability that the site x + Sn is unoccupied, that is, that 
U**{x + Sn) = 0, is equal to the probability that none of the branching ran- 
dom walks [/? succeeds in placing a particle at location x at time n. Since the 
attached branching random walks are independent of the random walk tra- 
jectory {Si}i<:n and the displacement random variables ^j, this probability 
is 

(6.2) P{site {Sn + x) vacant} = (1 - Ui{x + Si+i + ^j)) 

i 

where n„ is the hitting probability function 

(6.3) Unix) := P{Unix) > 1}. 

Proposition 31. There exists C > such that for all n > 1 and all 

sites X S Z^, 

(6.4) Unix) > — ^^"^ 



C + ^ log n 



Proof. By the fundamental identity, EUnix) = Pnix). By the second 
moment formula (5.7) of Proposition 25, 



n-1 



(6.5) EUnix)'' = Pnix) + J^P,(Z)P2_,(X - z) 

i=0 z 
n-1 

<Pnix) + Y,Yl Piiz)Pn-iix - z)Pn-^iO) 



1=0 z 

n-1 



Pnix) + Pnix) ^Pn-^iO) 



1=0 



< Pnix)iC + Alogn) 



imsart-aop ver. 2009/08/13 file: BRWpublish. tex date: April 8, 2010 



BRW OCCUPATION STATISTICS 



31 



Here we have used the fact (Lemma 11) that Pn-i{x) is maximal at the 
origin x = 0, together with a strong form of the local central limit theorem 
(specifically, the fact that the error in the local limit approximation is of order 
0(n~^), which is summable). The result (6.4) now follows immediately from 
the Cauchy-Schwartz inequality P{X > 0} > {EX)'^ / EX"^ , valid for any 
nonnegative random variable X. □ 

The lower bound (6.4) leads easily to a useful upper bound for the prob- 
ability that site x is vacant. Partition the indices i < n into two sets, the 
good and the bad indices, as follows: Fix a large constant k < oo, and say 
that index i is good if |S'i+i + < k\/1; say that i is bad otherwise. By the 
local central limit theorem, there is a constant C > not depending on k 
such that for every good index i > 

(6.6) Pi(x + 5i+i+e.) >C"e-2«V^. 

Thus, relations (6.4)-(6.2) and the concavity of the logarithm function imply 
that for a suitable constant C" > not depending on k, 



3.7) P{site {Sn + x) vacant} < exp < 



2 



i log i 

good, \x\^<i<n 



Lemma 32. Let {in} be any sequence of real numbers such i/iai lim„£„ = oo 
and lim„ log £„/ log n = 0. Then for every b > and every e > there exists 
K sufficiently large that 



lim sup P 



y 

^ i log i 




i good, £^<i<n 

Proof. The hypotheses regarding the growth of ensure that 

n 

Ln ■■= ^ l/{ilogi) — > oo. 

i=ll 

Hence, it suffices to show that for some < < 1, if «; is sufficiently large 
then 



(6.9) P\ V -TT^>gLr 

• ^-^ I log I 



< e 



i bad, £'}. <i<n 
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for all large n. Recall that an index i is bad if jiSi+i > kv^. Cliebyshev's 
inequality implies that for any e > 0, if k is sufficiently large then P{\Si^i + 
> Ky/i} < e^; hence, for large n, 

^ ■^{\S^+lH^\>'iVi} ^ ^3^ 

ll<i<n 



^ i log i ' 



It now follows by the Markov inequality that 
(6.10) P < 

The relations (6.10) clearly implies (6.9), and therefore prove (6.8). □ 

Proof of Proposition 30(A). For any e > 0, inequality (6.7) and 
Lemma 32 imply that for all large n, for any displacement x of magni- 
tude < In, the probability that site (x + S'n) is vacant is less than 1e. There- 
fore, the expected number of vacant sites in the ball B{Sn',in) given the 
event Gn is, for large n, no larger than Anei'^. The assertion (A) of Theo- 
rem 8 follows directly, by the Markov inequality. □ 

6.3. Proof of Proposition 30 (B). The second assertion (B) of Proposi- 
tion 30 can be proved in virtually the same manner as Theorem 6. Following 
is a brief sketch. Set 

Wn := ^ particles of generation n within distance of Sn 
in the size-biased BRW. 

By representation (6.1), 

n-l 

(6.12) Wn = 2 + Y, E Ul{x + Si+i + i^), 

i=l \x\<l^ 

where Uj{x), Sn, and satisfy conditions (Ta)-(Tc) of section § 5.3. The 
distribution of the sum on the right side is analyzed by decomposing it as 
r„ -|- A„, where now 

n 

(6.13) r„ ^ /^,(x + 5») and 

i=2 \x\<e„ 

^"•=E ( E iUizl{x + Si + C^.i)-Pix + S^)) \ . 
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By calculations similar to those used in proving Lemmas 27, one shows that 
(6.14) lim ET n/{Trillogn) = A/2; 

n— ^-oo 

lim Var(r„)/(7r^2 ^^g^) ^ Q; 

n— !>oo 

lim Var (A„)/(7r^^ log n) < and 



n— >oo 

EAr, = for all n > 1. 



Given these estimates, one now obtains the desired conclusion, that Wn is 
of order Op{£f^ logn), by the same simple argument as in section 5.6. □ 
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7. Occupied Sites in Dimension 2. 

7.1. Hitting probability function. For simplicity we consider in tliis sec- 
tion only the binary fission case; the case of a general offspring distribution 
with mean 1 and finite variance can be handled similarly. The proof of The- 
orem 7 will be based on careful analysis of the hitting probability function 
Un{x) defined by equation (6.3) above. The connection with the total num- 
ber f2„ of occupied sites at time n is obvious: EQn = X^a;^n(2;)- Thus, our 
goal will be to bound the function Un from above. (A good lower bound has 
already been obtained in Proposition 31.) Our main result is the following 
proposition. 

Proposition 33. There exist constants Ci,C2 < oo such that for all 
n> 2 and all sites x G Z^, 

Ci ( 

(7.1) Un{x) < — exp -C2 , 

nlogn \ n J 

and hence for some C > we have that 

C 



(7.2) E^n = 



log n 



Theorem 7 follows as a direct consequence of (7.2) and Kolmogorov's 
estimate (1.1). 

To obtain upper bounds on the function Un{x), we will exploit the fact that 
it satisfies a parabolic nonlinear partial difference equation. Recall that P is 
the Markov operator for the simple random walk, that is, for any bounded 
function w : Z'^ ^ R, 

z-x£j\f 

Lemma 34. Assume that the offspring distribution is double- or-nothing. 
Then for each n > and each x ^U^, 

(7.3) Un+l{x) = ^Un{x) - ^(Pn„(x))2. 

Proof. The event {Un+i{x) > 0} can only occur if the first generation 
is nonempty, and hence consists of two particles with locations in M. This 
happens with probability 1/2. One or both of these particles must then 
engender a descendant branching random walk that places a particle at site x 
in its nth. generation. Since the two descendant branching random walks are 
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independent, with starting points randomly chosen from A/", this happens 
with probabihty 2p{l — p) + p^, where p = Pti„(x). □ 

To extract information from the nonhnear difference equation (7.3) we 
will use the following standard comparison principle. (Compare, for example. 
Proposition 2.1 of Aronson and Weinberger (1975).) 

Lemma 35. Let Un{x) and Vn{x) he functions taking values between 
and 1 that satisfy the following conditions: 

(7.4) Un+i{x) = Fun{x) - ]^{Fun{x)f and 

(7.5) Vn+l{x) > FVnix) - ^iFVn{x)f. 

If Vo{x) > Uo{x) for all x, then 

(7.6) Vn{x) > Un{x) for alln>0 andx € 

Proof. Define A„(x) = Vn{x) — Un{x); then by the hypotheses (7.4)- 
(7.5), 

(7.7) An+i{x) > PA„(x) - -(Pn„(x) + ¥vn{x))FAn{x). 

Since ti„ and Vn take values between and 1, so does the average (Pn^ + 
Pt;„)/2. Therefore, (7.7) and the induction hypothesis imply 

A„,+i(x) > PA„(x) (^1 - ^ {Funix)+Fvnix))^ > 0. 

□ 

The trick is to find a function u„ that satisfies inequality (7.5) and domi- 
nates uq. To this end, fix k > and define 

(7.8) Vn{x) = — exp f- ^"J^^ 

n log n \ 2n 

where 

/3„ = /3(l--^) and /3 = 5/2. 

Lemma 36. There exist A'^o G N and kq independent of Nq such that 
for all n > kq and n > Nq, 

(7.9) Vn+i{x) > Fvnix) (l - ^Fv^ix) 
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The (rather technical) proof is deferred to section §7.3 below. (See Bramson et al. 
(1993) for a similar argument in the context of the KPP equation.) Given 
Lemma 36, Proposition 33 is an easy consequence. 

Corollary 37. There exist A''i G N and k> such that for all n>0, 

(7.10) Unix) < VNi+n{x) < 1. 

Proof. Choose A^i > Nq such that k := iVilogiVi > kq. For such a 
choice of (iVi , k) we have 

Uo{x) = l{a.=o} < vnA^) < 1- 

Moreover, by Lemma 36, the function Vn{x) := Vn+Niix) satisfies (7.9). The 
conclusion now follows from the Comparison Lemma 35. □ 

7.2. Representation of the conditional distribution. Revesz Revesz (1996) 
considers a branching random walk on that is identical to the branching 
random walk we have studied, except that the particle motion is by Gaussian 
A^(0, 1) increments rather than Uniform-AA increments. One of the main 
results of Revesz's article asserts that, conditional on the event that there is 
at least one particle of the nth generation in the ball B of radius g = vr"^/^ 
centered at the origin, the expected total number of such particles is of 
order G(logn). His argument seems to rest on the (unproven) assertion (see 
the first two sentences of his Proof of Theorem 3) that conditional on the 
event that a region C is occupied by at least one particle at time t, the 
branching random walk consists of a single pinned random walk off of which 
independent branching random walks are thrown. There is no proof of this 
assertion (in fact, it is not even stated clearly, as far as we can see). 

We believe that Revesz' assertion is false. The purpose of this section is 
to give a representation related to that of Revesz' for the conditional law of 
the occupation random variable Un{x) given the event 

:= {Unix) > 0}. 

This representation is similar to Revesz' in that it consists of independent 
branching random walks thrown off a random path from (0,0) to (n,x); 
however, the distribution of the random path is not that of a pinned simple 
random walk, but rather that of a u— transformed simple random walk. This 
is defined as follows: 

Definition 38. For each site x and integer n > 1 such that Un{x) > 0, 
the u— transformed simple random walk with endpoint {n,x) is the n— step, 
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time-inhomogeneous Markov chain {Xm}o<m<n on Z'^ with initial point 
and transition probabihties 

(7.11) qUz, y) := P{X^ = y \ X^^, = z) = Pi{y - z) """-"^"^ " 



¥Un-m{x - z) ' 



Remark 39. Except in the trivial case n = 1, a transformed random 
walk is not a Doob /i— process, because the hitting probability function Un{x) 
is not space-time harmonic for the simple random walk, by equation (7.3) 
above. But a pinned random walk is an /i— process: In particular, the one- 
step transition probabilities of a pinned random walk conditioned to end 
at are given by 

(7.12) C(^,y) = Pi(y-^)-^"-'"^'^""^^ 



Pn—m+l{,Xn 

Since the function Pn-m{z,Xn) is space-time harmonic, the transition prob- 
abilities q* are not the same as those of the n— transformed random walk. 

Lemma 40. If Un{x) > then the u— transformed simple random walk 
with endpoint (n, x) is well-defined, and with probability one ends at X^ = x. 

Proof. What must be shown is that the Markov chain with transition 
probabilities (7.11) will visit no states (m,z) at which the denominator 
¥un-mix — z) is zero. This is accomplished by noting that as long as X^^i is 
at a site z such that Un-m+i {x — z) > 0, then by Lemma 34 the denominator 
¥un-mix — z) > 0, and so there is at least one site y among the nearest 
neighbors of z such that Un-m,{x — y) > 0. By (7.11), the next state X^. will 
then be chosen from among the nearest neighbors such that Un-m{x — y) > 0. 
This proves that the Markov chain is well-defined. The path ends at Xn = x 
because is the only site at which uq > 0. □ 

Our representation of the conditional distribution of the random variable 
Un{x) given the event Gn^x requires four mutually independent sequences of 
random variables: 

{Ua) {Xra}Q<rn<n is a u— transformed simple random walk with endpoint 

{n,x); 

(Ub) {^m('W^)}o<m<n; weZ'* independent Bernoulli(/3m(ty)) random vari- 
ables; 

(Uc) {U^{y)}i>o are independent copies of the branching random walk {Um{y)}', 
and 

(Ud) {£,i}i>0 a-re independent and uniformly distributed on M. 
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The Bernoulli parameters are 



note that for large values of n — m the parameters Pm{w) are uniformly 
close to 1/2, because Un-m{x — w) is bounded by the probability that the 
branching random walk will survive for n — m generations. 

Proposition 41. Assume that the offspring distribution is double-or- 
nothing, and let x be a site for which Un{x) > 0. Then 



Proof. The assertion (7.14) is equivalent to the assertion (Claim 42 be- 
low) that the conditional distribution can be simulated by the following 
Method A: (1) Let a particle execute a n— transformed simple random 
walk {Xm}m<n with endpoint {n,x). (2) At each location {m,Xfn), where 
< m < n, toss a /3m(-^m) — coin to determine whether or not to attach a 
descendant branching random walk. (3) On the event that the coin toss is a 
Head, create a new particle let it make one jump ^m+i to a neighboring 
site, and then attach an independent branching random walk starting from 
this new location. (4) Count the total number of particles, including ^, that 
land at site x at time n. 

Claim 42. This simulates the conditional distribution of the total number 
of particles at site x in generation n given the event {Un{x) > 1}. 

This claim is proved by induction on n. The case n = 1 is routine, but for 
the reader's convenience we shall present the argument in detail. First, the 
only sites x such that ui{x) > are the nearest neighbors of the origin, so 
we assume that x is one of these five points. Since uq = 6q is the Kronecker 
delta function, Pno(x) = 1/5, and so /3o(0) = 1/(2 - 1/5) = 5/9. Now 
consider the first generation Zi of the branching random walk: this will be 
empty unless the initial particle fissions, in which case the two offspring are 
located at randomly chosen nearest neighbors of the origin. Consequently, 



(7.14) 
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the unconditional distribution of Ui{x) is 

m(^) = o}4 + ixlxi 

P{C/i(x) = l} = ix2xlxi; 
P{C/i(x) = 2} = ixixi 

It follows that the conditional distribution of Ui{x) given the event {Ui{x) > 
0} is that of 1 plus a Bernoulli(l/9) random variable. This coincides with 
the distribution of the random variable produced by Method A, because 
-Bo(O) = 1 with probability /3o(0) = 5/9, and on this event the particle 
jumps to X with probability 1/5, leaving a second particle at x. 

Next, consider the branching random walk conditioned to have at least one 
particle at site x in generation n > 2. The first generation must consist of 
two particles, at least one of which produces a descendant branching random 
walk that places particles at x in its (n — l)st generation. Conditional on the 
event that two particles are produced by the initial particle (that is, the event 
{Zi = 2}), each will have chance p := ¥un-i{x) of producing a descendant 
at site X in generation n; consequently, the conditional probability that both 
particles will do so, given that at least one does, is 

2 

^ p/3n-l(0). 



p2 + 2p(l-p) 

Moreover, given that either one of the particles produces a particle at site x 
in generation n, the conditional probability that its first jump is to site y £ M 
is 

f7-tK\ P f ^ Un-l{x-y) 

^'■''^ "^'^'^ Pn„_,(x) ' 

this is the distribution of the first step of a n— transformed random walk with 
endpoint {n,x). Thus, a version of the random variable Un{x), conditional 
on {Un{x) > 1}, can be produced by the following two-step procedure: 

(1) Place a particle 77 at a randomly chosen neighbor y of according to 
the distribution (7.15), and attach to it a branching random walk conditioned 
to produce at least one descendant at site x — y in its (n — l)st generation. 
By the induction hypothesis, the contribution of offspring of r] to site x in 
generation n will be 

(7.16) 1 + Sm,(^m)f/r-m-l(^ -Xrn- U+l)- 

m.=l 
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(2) With probability p/3„_.i(0), do the same with a second particle r. 
Observe that, conditional on the event that this second particle r is attached, 
the contribution to site x in the nth generation will have distribution 

C {Un-l{x - X,) I Un-l{x - X,) > 0) . 

Since the particle r is attached with probability p/3„_i(0), where p is the 
probability that a particle born at time will put a descendant at site x 
in generation n, step (2) has the same effect as this alternative: (2') With 
probability /3„_i(0), place a second particle r at a randomly chosen (that 
is, uniformly distributed) neighbor y of 0, and attach an independent copy 
of the branching random walk. This, together with the representation (7.16) 
of the number of offspring of r/ at x in generation n, shows that the total 
number of particles at x in generation n will be 

(7.17) 1 + Bm{X^)U::i^_,{x - X,n - U+l), 

as desired. This completes the induction argument, and thus proves (7.14). 

□ 



7.3. Proof of Lemma 36. Let x = (xi,X2); then 

5 

where 
Then 

1 



Vn+Xix) -^Vnix) + U^Vn{x)f 



^ ^ 5 V (n + 1) log(n + 1) ^ V 2 y ^ 



10 

where 

r, /3n /3n+l 
Pn — • 

n n + 1 



g-/3n/{2n) \ 

H 77; Vn{x)wn{x) 
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Therefore it suffices to show that there exist Nq and kq independent of Nq 
such that for all n > kq a-^d for all n > Nq, the following holds: 

(7.18) (n + l)log(n + l)^"PV 2 J ""^"^ 



e 



-/3n/(2n) 



+ Vnix)Wnixf > 0. 



a. Estimate of e^-A^") : First, 

nlogn ^ (n + 1) log(n + 1) — nlogn 



(n+l)log(n + l) 



(n 


+ 1) log(n + 1) 


1 


nlog(l + 1/n) 


n + 1 


(n+l)log(n + l) 


1 


1 


n + 1 


(n + 1) log(n + 1) 



1 



1 — - . . ... T+O 



1 + ^ + — 

n n log n 



1 



Therefore, recall that /3„ = /3(1 — 1/logn), 



p/3„/(2n)_ 



n log n 



(n + l)log(n + 1) 

zn inlogn \n-^ J J \ n nlogn 

=1 + _ /? + 2 ^ /3g-4/3 + 8 ^ ^ M 

2n 2nlogn 8re^ 

Since /3„ — )• /3 = 5/2, there exists Nq G N such that for all n > Nq, 

(7.19) e/^-A^"), >l+^^-^+A>l- 

(n + 1) log(n + 1) 2n 2?ilogn 8n^ 

b. Estimate of Since /3„ = /3(1 — 1/logn), we have 

^ _^ / I -1/logn l-l/log(n + l) 
" \ n n + 1 

^1 — (n + l)/logn + n/log(n + 1) 
n{n + 1) 
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However, 

n+1 n log(n + 1) + n log(l + 1/n) 



logn log(n + l) lognlog(n+l) 



logn ylognlog(n + 1) 
so it follows that 

(7.20) On = ^ 1 - V ^ + 0( VQog ^ log(^^ + 

" n(n + 1) 

Claim 43. Enlarging Nq if necessary, we have that for all n > Nq, 

0^ 1 

(7.21) f39n-^> 



log n 

Proof of the claim. Since /3„ = /3(1 - 1/logn), 



n+lj \ log 77, ylog 77 log(77 + 1) 



2 1 
+ 



log 77 (log 77) 



^ log 77 \ log 77 ^ 

The relation (7.21) follows since /? = 5/2 > 1. □ 

c. Proof of (7.18) for \x\ > 3n: (7.21) implies, enlarging Nq if necessary, 
that for all 77 > A'^o, On > 2/77^. Hence when \x\ > 377, 

0n\x\^/2> Pn\Xi\/n, 7=1,2, 

and 

5exp ( — - — I > Wn{x). 
The relation (7.18) follows by noting (7.19). 
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d. Estimate of Wn{x): For \x\/n sufficiently small, Taylor expansion yields 



2n 2n log n 8n^ 



+ 2 + 



+ 



(7.22) 



12n4 

4^4 



+ 2 + 



+ 



12n4 



+ r^] 



n 



X2y 



=5 + 



+ 



/3 



^ 

2n 2n log n 



+ 



/32 
8n2 



n / 



+ 



+ 



/34(x| + x|) 



+ 



e. Estimate of e^"/^^") 
n>iVo, 



12n4 

n log rt r^^r-r^(0 1^1 2 



„/3n/{2n) 



(n+i)log(n+i) • exp(0„,|x|V2): By (7.19), for all 
n log n 



(7.23) 



> 1 + 



(n + l)log(n + l) 
13-2 13 + 2 2 
2n 2n log n Sn^ 



• exp 



2 

1 + 



+ 



>1 + 



/3 + 2 



2n 2?i log n 
f: Their difference: By (7.23) and (7.22) 



+ 



el\x\' 



5e 



>5 



/3„/{2n) . nlogn 



{n + 1) log(n + 1) 



exp 



Wn{x) 



(7.24) 



+ 



13-2 _ (3 + 2 
2n 2n log n 

8n2 8n2 



/5 



^ 

2n 2n log n 



+ 



2U|4 



Since /3 = 5/2, 
(7.25) 5 



+ 



/3 



59i\x 



x\ 
12n^ 



_ (3 + 2 
2n 2n log n 



/3 



l_ 

2n 2n log n 



10 



n log n ' 
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and enlarging A'^o if necessary we can assume that for all n > Nq, 

Moreover, 9n ~ it follows that for all n sufficiently large, 

^'■^'> 8 12n4 - I 8 12n4 j "l"^! > 2n4 ' 



g. Proof of (7.18) for Sn > \x\ > vTOn, where 5 > is sufficiently 



small: By (7.21), when |x| > VW 



v? J n log n 

Hence, by (7.24), (7.25), (7.26) and (7.27), the relation (7.18) holds for x 
such that \x\ > i/TOn and \x\/n sufficiently small. 

h. Proof of (7.18) for 3n > \x\ > 5n: By (7.21), for all n > Nq, 

/32 1 
~ /3n^ log n 

Hence when |x| < 3n, 



exp I — — '— I > " \ > exp ( — I • exp 



x\ 



expi^^^) V logn; V 5n2 

^ ori'^ log n y 

By (7.19), to show (7.18) it is sufficient to show that for all n sufficiently 
large, 

g/3„/(2n) _|_ g-/3„a;i/n _|_ ^P„xi/n _|_ ^-l3„X2/n _|_ ^l3„X2/n 



< 5 exp — • exp 



log n J \ 5?2^ 

Since < 3n, 



1 - exp • exp , 

log n / / \ on^ 
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Hence it suffices to show that 



hminf inf < 5exp , 

(7 28) " 3n>|a;|><5n [ \ 5re^ 

_ _j_ Q-l3nXi/n _|_ ^PnXi/n _|_ ^-l3„X2/n _|_ g/3„a:2/n^ | > 0. 

By elementary calculus, 

g-/3na:i/n _|_ g/3„a:i/n _|_ g-/3„a:2/n _j_ g^„a:2/n < 2 + g"''"!^!/" + g'^"!^!/", 

thus 



5n^ 

>5exp (^1^) - 3 - e-^"l-l/" - e/^-l^l/". 

Relation (7.28) now follows from the simple fact that 

f{x) := 5e^'/^ - 3 - - e"^ 

is strictly increasing for x > and equals only when x = 0. 

i. Proof of (7.18) when |.x| < VTOn: Since \x\ < \/TOn = o(n), by relations 
(7.24), (7.25), (7.26), (7.21) and (7.27), we need only show that there exists 
kq such that if k > kq and n > Nq, then 



e 



-/3n/(2n) 10 



Vnix)Wnix) > 



10 n log n 

Since Wn{x) > 5, when |x| < \/10re, 

— Vn{x)wn{x) > -— j exp(-6/3), 

lU lUnlogn 

so Ko can be chosen as 4exp(6/3), which is independent of Nq. 



□ 
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